Stability Spectroscopy of Rotons in a Dipolar Bose Gas 
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We study the stability of a quasi-one-dimensional dipolar Bose-Einstein condensate (dBEC) that 
is perturbed by a weak lattice potential along its axis. Our numerical simulations demonstrate that 
systems exhibiting a roton-maxon structure destabilize readily when the lattice wavelength equals 
either half the roton wavelength or a low roton subharmonic. We apply perturbation theory to 
the Gross-Pitaevskii and Bogoliubov de Gennes equations to illustrate the mechanisms behind the 
instability threshhold. The features of our stability diagram may be used as a direct measurement 
of the roton wavelength for quasi-one-dimensional geometries. 



It is widely believed that ultracold, gaseous samples of 
bosonic atoms or molecules possessing sufficiently large 
dipole moments will exhibit internal structure reminis- 
cent of the roton in superfluid helium [l[ . The basic phe- 
nomenology of the roton, in analogy with helium, is a 
local minimum in the quasiparticle dispersion uj(k). The 
existence of such a minimum is predicted to lead to a 
host of attendant phenomena in these dilute gases, in- 
cluding structured ground-state density profiles 4L re- 
duced and anisotropic critical superfluid velocity [1 
enhanced sensitivity to external perturbations [7[ , abrupt 
transitions in Faraday patterns [1, [j| , short- wavelength 
immiscibility phases [lOj, and strongly-oscillatory two- 
body correlations on the roton length scale Signa- 
tures of the roton in Bragg spectroscopy of trapped dipo- 
lar Bose-Einstein condensates (dBECs) have been calcu- 
lated in Ref . [l2| • While these many exciting predictions 
are in principle observable in current experiments with 
highly magnetic atoms (T34l5| . or in future experiments 
with electrically polar molecules [l6l - [l8j . rotons remain 
to be seen directly in dBECs [l9j . 

In a dipolar condensate, the roton represents a mode 
of finite wavelength that has an anomolously low energy 
- hence the minimum in cj(k). The location of this min- 
imum is given by a momentum approximately equal to 
H/lt, where l t = y fi/mwt denotes the harmonic oscillator 
length of the tightest confinement of the dBEC (usually 
along the polarization axis). Without this confinement, 
a homogeneous dBEC would be energetically unstable 
to collapse due to the attraction between dipoles that 
are aligned head-to-tail. In the presence of this confine- 
ment, the collapse is prevented by the zero-point energy 
introduced by the confinement, at least up until a crit- 
ical dipole moment or density. When this critical pa- 
rameter is exceeded, the condensate collapses in localized 
"clumps" of size A rot ~ It, that is, via a dynamical insta- 
bility into the roton mode [2(j ■ For dBECs that are just 
barely stable, the energy minimum w(fc rot ) is achieved at 
the roton wave vector fc rot = 27r/A ro t. 

A low-energy roton is therefore a mode that is linked 
intrinsically to condensate instability. This link sug- 
gests that the dBEC may respond nontrivially as an ob- 
ject of conventional spectroscopy, that is, that it would 
absorb and distribute energy differently from different 



wavelengths of an applied probe. For example, if one im- 
poses on a dBEC a weak potential of periodicity Al (via 
a ID lattice beam [2l|, for example), then one expects 
to trigger an instability most easily when A L « A rot . In 
this Letter we verify this conjecture via mean-field sim- 
ulations of a quasi-lD (qlD) dBEC. Strikingly, we find 
structure in addition to this main peak, analogous to 
nonlinear spectroscopy of atoms or crystals in strong laser 
fields. Namely, we see hastened destabilization by probes 
with wavelengths that are integer multiples of A ro t , rem- 
iniscent of "multiphoton" scattering. We also observe a 
similar feature at the shorter wavelength Al = A ro t/2, 
reminiscent of resonant Raman coupling. Since the sta- 
bility of a condensate is easy to assess experimentally, 
the observation of such structures constitutes a direct 
sta&iZi£j/-spcctroscopic measurement of the roton wave- 
length as well as a novel signature of roton physics. 

Consider a dBEC that is tightly confined in the y and 
z directions by a harmonic trap of frequency u>t, with 
no trapping potential in the x direction. The dipole mo- 
ments of the constituent atoms or molecules are polarized 
along z. To this initially-stable system is applied a probe 
in the form of an optical lattice potential 



U{x)/tvjj t 



(1) 



whose periodicity Al = 27t//cl and (dimcnsionlcss) lattice 
depth s = —Rc{a(uj)}Io/2eochujt > are tunable pa- 
rameters. Such a potential can be realized using retrore- 
flected or crossed off-resonant laser beams of peak in- 
tensity Iq [2l[. We take the single-mode approximation, 
assuming that the order parameter is a Gaussian of width 
it = y/ h/mu>t in the directions of tight confinement [22| . 
We suppose that there are N = 2Lnm atoms spread 
over the periodic domain x £ [—L,L], where nuj is the 
one-dimensional integrated density. Expressed in natural 
length (l t ) and energy (fkut) units, the Gross-Pitaevskii 
equation for this situation reads 



fitp(x) = - -d 2 x i'{x) + U(x)i/)(x) 

+ N [ dx'il)*(x')ip(x')V(x - x')ip(x). 



(2) 



FIG. 1: Stability spectrum for a qlD dBEC in a weak lattice. 
The experimental parameters are for 164 Dy trapped at u)t = 
2ir x 5 kHz with density uzd ~ 10 15 cm -3 and dipole moment 
d — 10/j.b- In Case I (II), the scattering length is tuned to 
a s — — 22ao (a s = — 24.5ao), yielding a roton wavelength of 
900 nm (850 nm). The vertical dotted lines denote A ro t/2, 
A ro t, and 2 Arot, which identify the spectral features. 



The momentum-space form of the qlD interaction po- 
tential is 



V{k) =2a s + 4a dd 

CO 
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(3) 



where g(q) = qcxp(q 2 )crk(q). Contact interactions con- 
tribute to V via the scattering length a s , and dipolc- 
dipole interactions contribute via the dipole length add = 
md 2 /3h 2 , where d is the dipole moment of the bosons. 
We remark that V(k) decreases monotonically as a func- 
tion of k. The qlD chemical potential \i is related to 
the three-dimensional chemical potential \xzd via fi = 
fj,3D — 1. In the absence of the probing potential U (x), the 
chemical potential and order parameter are easily found 
to be fi^ = 71idV(0) and tp^ = 1/V2L, respectively. 

To determine the stability of the dBEC for a given 
lattice periodicity Al and intensity s, we solve Eq. 
and determine the excitation spectrum within the Bo- 
goliubov approximation. Dynamical instability is indi- 
cated by a complex frequency in a given mode, which 
also identifies the mode that triggers the instability. Nu- 
merical results are plotted in Fig. [TJ which shows two 
sample stability spectra for dBECs consisting of 164 Dy 
atoms. Plotted is the intensity s versus periodicity Al 
of the probe. The shaded region below each curve repre- 
sents the stable region, since higher intensitites perturb 
the condensate more strongly and lead to collapse. In 
one case (I, yellow region) the scattering length is as- 
sumed to be a s — — 22ao, while in the other (II, brown), 
it is a s — — 24.5ao- These differences correspond to dif- 



FIG. 2: (Color online) The green (dotted) and blue (solid) 
curves are respectively the Hartree-Fock and Bogoliubov spec- 
tra for the particular cases considered in Fig. [T] 



ferent roton spectra, as shown by the blue lines in Fig. 
[5] Case II has a "softer" roton than case I, i.e., a lower 
excitation energy w(fc ro t)- Along the lines of Ref. [23|, we 
have checked that our quantum and thermal depictions 
for both case I and case II, at an assumed temperature 
of 50 nK, remain small compared to the condensate frac- 
tion. Similar results can be found for all other species of 
strongly dipolar condensates, although there is a trade- 
off between density nnj and dipole length add for a given 
system to rotonizc. 

In both cases, the spectra clearly show three features, 
which arc more prominent for softer rotons: i) A weak 
laser drives the condensate to instability for probe pe- 
riodicities near the roton wavelength, accounting for the 
lowest dip in Fig. [TJ Indeed, the weakest laser required to 
de-stabilize the condensate has periodicity Al = A ro t . We 
emphasize that at the threshhold of stability, our simula- 
tions show that the lattice probe is not chopping the con- 
densate into smaller, unstable pieces, but rather is weakly 
perturbing it. ii) A secondary dip in stability appears 
when the lattice spacing is twice the roton wavelength. 
In terms of the analogy to optical absorption in a nonlin- 
ear medium, this situation is reminiscent of the excita- 
tion of the resonant state by absorption of two photons of 
half the energy required to excite the state directly, iii) 
Finally, a third dip in stability is observed at half the ro- 
ton wavelength. We have also observed these structures 
in full 3D simulations of horizontal, cigar-shaped dBECs 
where the roton "wavelength" is less well defined. 

The origin of these features becomes clear within per- 
turbation theory. Rigorous perturbation theory for the 
Gross-Pitaevskii equation was worked out in Refs. [Timl- 
in . Let ip(°) and represent the order parameter and 
chemical potential in the absence of the probe U(x). Ex- 
panding ijj = -0'°^ + ipW, f-i = /J*- ' + /U^ 1 ', one derives the 
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first-order perturbation equation 

^V 0) = -\d 2 x i^ + U4^+2n 1D J dx'^V(x-x'). 

(4) 

Given the condition (rj)W = 0, the solutions to (U]) 

are easily found to be = s/2 and 



in 



COs(fcL^)- 



(5) 



We have introduced the Hartree-Fock energy [27[ 
e(k) = k 2 /2 + 2niDV(k). Figure [2] compares e(k) to the 
Bogoliubov spectra uj(k) = y/k 2 /2 (fc 2 /2 + 2n 1D V(k)) 
for the two cases of dBEC considered in Fig. Q] 

The Hartree-Fock spectrum of a soft-roton dBEC 
posesses a shallow minimum at k « fc rot , and this min- 
imal energy approaches zero as the roton softens. The 
effect of this low-lying mode is to strongly perturb the 
ground state, as suggested by JS]). Indeed, a modest 
perturbation oc s becomes amplified in the density by 
a factor s/e(fc ro t)- The m th order of perturbation the- 
ory additionally introduces density modulations of the 
—r—j — v cosimkj.x). That is, the 



general form n\^{x) 



roton can be driven by overtones of the fundamental wave 
number of the lattice, when mk^ ~ A: rot , or Al ~ mA ro t. 
The stability minimum corresponding to m = 2 can be 
seen in Fig. [T] Even away from one of these resonant 
conditions, the probe laser introduces density modula- 
tions that manifest themselves in a mean-field potential 
U m {(x) — N J dx'\ip\ 2 V(x — x 1 ) which combines with the 
probe field itself to make a combined potential 



U, 



(x) =U(x) - fi + U m f(x) 



4e(fc L ) 



cos 



<fc L x) + 0{s 2 



(6) 



Thus again, whatever influence the periodic probe poten- 
tial has on propagation of the excited states in the con- 
densate, this effect is amplified for probes in the vicinity 
of the roton (or a subharmonic thereof, as a higher-order 
effect). 

There remains the task of explaining the feature in 
the stability spectrum at Al ~ A ro t/2. This cannot be 
done by considering the ground state only. Rather, it is 
necessary to explore how the spectrum of excited states is 
modified by the combined potential energy U c ©. The 
Bogoliubov de Gennes (BdG) equations describing the 
excitations can be written compactly as 



H 



-X 



X 



-H 



X 

fi-C-X 



E 



(7) 

where Ho = —\d 2 + U(x) is the free-particle 
Hamiltonian, C[x](x) = U m [(x)x{x) describes di- 
rect interactions with the mean-field, and X[x](x) — 
N J dx'x(x')ip(x')V(x — x')ip(x) is an integral opera- 
tor describing exchange interactions. The functions 



Uj(x) and Vj(x) define the usual Bogoliubov trans- 
formation of the quantum fluctuation field operator 

5$ = V, ajUj + a]v 



and they allow one to write 
the nontrivial part of the grand canonical Hamiltonian in 
approximate diagonal form as H — fiN ~ ■ Ej&j&j. In 
the absence of the perturbation, Eq. ([7]) is easily solved 
using complex exponentials parameterized by momen- 
tum fc, yielding the unperturbed energies E^ = uj(k). 
The perturbation ([!} separates modes of definite par- 
ity, so we write the unperturbed modes in a basis of 



cosine and sine functions: 



1 k,s 



(x) 



Uk 



sin(fcn)/ VI, 



u 

,(0), 



(0) 



(x) 



Uk cos 



{kx)/y/L, 



u k,'s( x ) = VkSm(kx)/VT 
are defined by Uk — 



J k ' c {x) = Vk cos(kx)/^/L, and 
The amplitudes u k and Vk 



k 2 /2+n 1D V(k) 
2w(fc) 



^ and Vk 



sgn(F(fc)) v /^Zl^mrT ; for fc > o ||. A system 

destabilizes when at least one of the excitation energies 
Ej vanishes. 

Assuming an initially stable condensate, the influence 
of U(x) on the excited states can also be approximated in 
perturbation theory. A perturbation theory for the BdG 
equations was developed rigorously in Refs. [HI, H(| [29| 
in the phase-density representation. Since we are chiefly 
concerned with first-order mode destabilization, we de- 
velop a tractable perturbation theory that paints a phys- 
ical picture of mode softening and naturally encompasses 
the degeneracy of the roton spectrum. After expanding 
Ek, Uk(x), and Vk(x) in perturbation series, and then 
substituting into Eq J7]), we derive the first-order per- 
turbation equation for the corrections E^\ u^\x), and 

fu(k) + \d 2 x - A<°) 

(uP+XW A« 
^ XV UP + A(D 

where X^[ X ] = Nj dx'x(x')^ m V{x - x') is the 

zeroth-order exchange integral operator, Uc (x) is the 
first-order combined potential given by Eq. ([6]), and 
A«[ X ] = N Jdx'x(x')^ 0) V(x-x') (^{x') +^ 1 \x)) 
is the first-order exchange integral operator. We can iso- 
late the energy shift E^' on the right-hand side by tak- 
ing advantage of the fact that the left-hand side vanishes 

whenever it is acted on by the operator J dx (uffi , v^j 

for any w(fc') = u(k). Degenerate perturbation theory is, 
in situations exhibiting a roton-maxon excitation spec- 
trum, generally necessary for a complete understanding 
of first-order effects; however, we will see that the stabil- 
ity structure at Al = A ro t/2 may be explained without 
accounting for degeneracy. To simplify our notation, we 
denote the hermitian matrix on the right-hand side of J5]) 
by A. 
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FIG. 3: Depiction of the overlap between the combined po- 
tential Uc(x) and the cosine and sine Bogoliubov modes for 
k = k L /2. 

It turns out that the matrix elements of A vanish in 
most instances, simplifying our analysis. Firstly, the 
sine modes completely decouple from the cosine modes 
as anticipated. Moreover, the matrix element of A be- 
tween any modes (uj^ , ujj.?^ and (u^\v^J vanishes 

unless the mode-matching condition \k' ± k\ = kt is sat- 
isfied. This follows from reasons of orthogonality, since 
all matrix elements are evaluated by integrating prod- 
ucts of three sine or cosine functions. Satisfying both 
the mode-matching condition and the degeneracy condi- 
tion w(fc') = uj(k) severely limits the number of nonzero 
matrix elements determining the modes (or degenerate 
mode mixtures) that shift at first order. 

The significance of the mode matching is illustrated 
schematically in Fig. [3] For any given perturbation 
probe with wave vector kt,, a standing matter wave of 
this wavelength is established, defining the combined po- 
tential U c (x) (upper panel). Against this backdrop are 
shown the cosine and sine modes of wave vector kt,/2 
(during some phase of their oscillation) in the lower two 
panels. We depict these modes because they always sat- 
isfy both the mode-matching and degeneracy conditions. 
The density oscillations of the cosine mode accumulate 
atoms at the maxima of U c (x), causing an increase in the 
mode energy by an amount E^\ 2 c = A^\ 2 k L /2- Con- 
versely, the density oscillations of the sine mode collect 
atoms at the minima of U c (x), thereby lowering the mode 
energy by the amount Ej^ , 2 s = —E^\ 2 c . Linear combi- 
nations of degenerate modes satisfying \k' ± k\ = kt, may 



also accomplish this spatial appropriation of atoms in an 
average sense, as the direct terms in the diagonal ele- 
ments appear as Ajj ~ / dx (^u^ 2 (x) + v^ 2 (x)j U c (x), 

where and Vj ^ are linear combinations of matched 
modes. This is the mechanism behind the first-order soft- 
ening of Bogoliubov modes and degenerate mode mix- 
tures, and it is similar to the discussion of "staggered 
modes" found in Ref. (30| . 

Given that the roton mode is a local minimum of the 
system dispersion, we expect to observe a lower stabil- 
ity boundary for all lattice periodicities in which the ro- 
ton softens to first order. This occurs when the roton 
mode- matches to itself (kt = 2fc rot ) and when it mode- 
matches to a degenerate phonon (kt fcrot)- These 
two scenarios correspond to the stability dips observed 
at Al = A ro t/2 and Al ~ A ro t, respectively. Of course, 
the latter structure was already expected due to strongly- 
enhanced, static density modulation, as discussed previ- 
ously. The structure at Al = A rot /2, however, occurs in 
spite of relatively weak density modulation. In that case, 
the roton responds directly to the perturbation, rather 
than to a strongly-amplified mean field potential. 

We note that our discussions of density response and 
mode softening remain applicable despite the increase 
in quantum and thermal fluctuations for rotonized sys- 
tems. There exists a minimum roton energy A below 
which either thermal or quantum fluctuations cause a 
strong reduction in the condensate fraction [23j . For un- 
perturbed dBECs with roton energies exceeding A, the 
perturbing lattice will likely drive the system towards 
an uncondensed or thermal state before reaching dy- 
namic instability. The criterion for the existence of a 
stable dBEC would then change from min{cj(k)} > to 
min{uj(k)} > A, effectively lowering the boundaries in 
Fig. Q] without changing the essential structure thereof. 

In summary, we have explored the novel stability fea- 
tures of quasi-one-dimensional dipolar condensates in a 
weak lattice. Rotonized systems destabilize rapidly when 
the lattice periodicity Al matches either half of the roton 
wavelength or a subharmonic of the roton. The observa- 
tion of these stability dips is a clear signal of rotonization, 
as well as a direct measurement the roton wavelength. 
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